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ABSTRACT 

We present one possible mechanism for the giant flares of the Soft Gamma-Ray Re- 
peaters (SGRs) within the framework of magnetar, i.e., superstrongly magnetized 
neutron star model, motivated by the positive period increase associated with the Au- 
gust 27 event from SGR 1900-1-14. From the second-order perturbation analysis of the 
equilibrium of the magnetic polytrope, we find that there exist different equilibrium 
states separated by the energy of the giant flares and the shift in the moment of inertia 
to cause the period increase. This suggests that, if we assume that the global recon- 
figuration of the internal magnetic field oi H > 10^^ G suddenly occurs, the positive 
period increase APt/Pt ^ 10^* as well as the energy ^ 10'*'* ergs of the giant flares 
may be explained. The moment of inertia can increase with a release of energy, because 
the star shape deformed by the magnetic field can be prolate rather than oblate. In 
this mechanism, since the oscillation of the neutron star will be excited, a pulsation 
of ~ ms period in the burst profile and an emission of the gravitational waves are 
expected. The gravitational waves could be detected by the planned interferometers 
such as LIGO, VIRGO and LCGT. 
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^ ■ 1 INTRODUCTION 



1.1 Giant flares of SGRs 



lO^*' erg s"* that emit brief 
10^^-10*^ erg s~*) recurrent bursts of soft (~ 30 keV) gamma-ray (see Hurley 1999 for a review). There 



Soft Gamma-ray Repeaters (SGRs) are persistent X-ray sources with X-ray luminosity ~ 10' 

(~ 0.1 s), intense ( 

are four known SGRs, SGR 1900-H4 (Mazets, Golenetskii & Guryan 1979), SGR 1806-20 (Laros et al. 1986), SGR 1627-41 

(Woods et al. 1999a) in the Galactic plane and SGR 0525-66 (Golenetskii et al. 1979; Mazets et al. 1979; Cline et al. 1980) in 

the Large Magellanic Cloud, and one possible candidate, SGR 1801-23 (Cline et al. 2000) towards the Galactic center. Their 

apparent associations with supernova remnants support that SGRs are neutron stars (Hurley 1999) . 

The spin periods and the spin-down rates of SGR 1900-1-14 (Woods et al. 1999b) and SGR 1806-20 (Kouveliotou et al. 1998) 
are measured from the quiescent X-ray pulses. If the spin-down is driven by magnetic dipole radiation, the implied magnetic 
fields are H J> 10*'' G, which are greater than the quantum critical field Her = m\c? /eh ~ 4.4 x 10*^ G (Baring & Harding 
1998; but see also Camilo et al. 2000). Although the spin-down may not be driven by the dipole radiation but a relativistic 
particle wind (Thompson & Blaes 1998; Kouveliotou et al. 1999; Harding, Contopoulos & Kazanas 1999; Rothschild, Marsden 
& Lingenfelter 1999), the rotational energy loss rate (|Tr2f2| ~ 10"^'* erg s~*) is not sufficient for the persistent X-ray emission, 
the particle winds and the soft gamma-ray bursts (^ 10^® erg s~*), and hence the superstrong magnetic field H J> 10** G 
is indispensable from energetics if the energy source is the magnetic field. Such superstrongly magnetized neutron stars are 
called "magnetars" (Duncan & Thompson 1992; Thompson & Duncan 1993, 1995, 1996), and SGRs are better understood 
within the framework of magnetar than other possibilities, such as accretion models (e.g., Kaplan et al. 2001; Liang 1995; 
Marsden et al. 1999; Cheng & Dai 1998, 1999). Anomalous X-ray pulsars (AXPs) are another group of X-ray sources that are 
similar to the SGRs but have no burst activity (see Mereghetti 1999 for a review). 

During the past ~ 20 years, two giant flares have been recorded from two of the SGRs: 1979 March 5 event from SGR 



© 2001 RAS 



2 K. loka 

0526-66 (Golenetskii et al. 1979; Mazets et al. 1979; Cline et al. 1980) and 1998 August 27 event from SGR 1900+14 (Hurley 
et al. 1999; Feroci et al. 1999; Mazets et al. 1999). These flares differ from the more common bursts in larger energies (<; 10*'* 
ergs), longer durations (^ 100 s) and harder initial spectra (~ MeV). In the August 27 event, the giant flare was intense 
enough to ionize the Earth's night-side lower ionosphere to the levels usually found only during daytimes (Inan et al. 1999). 
A radio afterglow was also found (Frail, Kulkarni & Bloom 1999). In the magnetar model, it is considered that the giant 
flares are produced by the global reconfiguration of the internal magnetic field, while the smaller bursts are produced by 
"crustquakes" in the neutron star (Thompson & Duncan 1995). 

Remarkably, during an 80 day interval in the summer of 1998 which contains the August 27 event, the average spin-down 
rate of SGR 1900-1-14 increased by a factor ~ 2.3 (Woods et al. 1999b). The samphng of the period history of SGR 1900-1-14 
is insufficient to distinguish between a long-term (~ 80 days) increase of the spin-down rate and a sudden increase in the 
spin period (a "braking glitch") at the time of the August 27 event .n In either case, this variable spin-down of SGR 1900-1-14 
should provide an important clue to the nature of the SGRs. 

Thompson et al. (2000) considered the physical mechanisms for the positive period increase of the observed magnitude 
APt/Pt ^ 10^^ directly associated with the August 27 event. They focused on two mechanisms: one is a loss of angular 
momentum by a particle wind and the other is an exchange of angular momentum between the crustal neutron superfiuid and 
the rest of the neutron star. The former mechanism can induce the observed spin-down only if an additional outflow of ~ 10*'' 
ergs Icists longer than the observed duration ~ 10* s and shorter than ^ 10^ s, and hence a fine tuning may be needed. The 
latter mechanism may cause the observed spin-down. However, at present, the exact cause for the positive period increase of 
SGR 1900-1-14 and the giant fiares remains unknown. 

In this paper, we will present other mechanism for the giant flare within the framework of magnetar, in which the global 
reconfiguration of the internal magnetic field of _ff ^ 10^^ G causes the giant flare of ^ 10** ergs as well as the positive period 
increase of APt/Pt ~ 10"*. 

1.2 Order-of-magnitude estimate 

Before discussing the details, we shall make an order-of-magnitude estimate. Let us consider a rotating star. There is a relation, 

j = in, (1) 

between the angular momentum J, the angular velocity Q = 2Ti/Pt and the moment of inertia T of the star. There are two 
ways to cause the spin-down: one way is to decrease the angular momentum and the other is to increase the moment of inertia. 
The mechanism considered by Thompson et al. (2000) is the former one. Our mechanism for the spin-down is the latter one. 
If the angular momentum J is conserved, the increase of the moment of inertia required for the positive period increase of 
SGR 19004-14 is of order AT/J ~ APt/Pt ~ 10"*. This corresponds to a rotational energy loss for the star \IQ.AQ,\ ~ 10*^ 
ergs, which is much smaller than the energy of the giant flare. 

Let us consider the internal magnetic fleld oi H ^ 10^^ G. This is somewhat larger than the external dipole fields H ~ 10^* 
G deduced from the period and the period derivative. Such larger internal fleld H ^ lO^'^ G is plausible if the magnetic fleld 
is generated by the helical dynamo in a new born neutron star (Duncan & Thompson 1992; Thompson & Duncan 1993). The 
magnetic field is a source of non-hydrostatic stress in the interior of the star. Hence the internal magnetic fleld induces a 
deformation of the star, and the magnetic deformation dominates the rotational one when H ^ 10'*(Pt/ls)"^ G. The induced 
difference in the moment of inertia of the star is given by order of 

AT (HV8^)(47rflV3) 4/ H \^/ fl \*/ M V' 

T ~ ~ GM^/R ~ V2xlO«G/ {w^cmj \IAMq J ' ^' 

where S is the ratio of the magnetic energy M ~ (ii'^/87r)(47r7?^/3) ~ 7 x 10*^(f//2 x 10^''G)^ergs to the gravitational energy 
|W| ~ GM^ / R ~ 5 X 10^^ ergs, and M and R are the mass and the radius of the star respectively (see below; see also Bocquet 
et al. 1995; Bonazzola & Gourgoulhon 1996; Konno, Obata & Kojima 1999). The deformation of the star also causes the 
gravitational potential energy W ~ —GM^/R to change by order of (Chandrasekhar 1969; Shapiro & Teukolsky 1983) 



*'l"'l/=-'°"(lir^)'(Kil;)'(dls)' -■^- P' 



AW " " ^*' « V^ « 

* Rothschild, Marsden & Lingenfelter (1999) argue that a sudden period increase at the August 27 event appears to be at odds based 
on the relatively large period derivative Pt ~ 1 X 10"^''s/s measured about three months before the giant flare (Kouveliotou et al. 1999) 
and just after the event (Murakami et al. 1999). However a reanalysis of all RXTE observation gives lower period derivatives, although 
the same data are used (Woods et al. 1999). Moreover, the timing noise of SGRs may be large (Woods et al. 2000) so that the long 
term trend may be more important than the short term trend. Therefore, we think that a sudden period increase is consistent with the 
observation. 
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Here we assume that the deformation is elhptical like a rotating star. Because the magnetic pressure is strongly anisotropic, 
this is a good approximation for globally coherent fields, such as a dipole poloidal field. In this case, 5 has the same meaning 
as the oblateness parameter e in equations (10.11.2) and (10.11.7) of Shapiro & Teukolsky (1983). The energy shift is of order 
S^ because the gravitational energy as a function of the oblateness parameter 5 is minimum aX 5 — Qx\ We will confirm that 
these crude estimates in equations (El) and (H) are good approximations in Section 4.2. 

Therefore, if the global rearrangement of the internal magnetic field H ~ 10^^ G occurs, the moment of inertia and the 
energy will change by order of AT/T ~ 10~* and AW ~ 10*^ ergs respectively, which is comparable with the required values 
for the positive period increase and the giant fiare of SGR 1900+14. In this way, the global reconfiguration of the internal 
magnetic field H ^ 10^ G may explain the positive period increase t^PtjPt ~ 10~ as well as the energy ^ 10 ergs of 
the giant fiare. This coincidence of the order-of-magnitude estimate and the observation may be by chance but is interesting 
enough to explore. Note that the energy source of this mechanism is not the magnetic energy but the gravitational energy. 
Furthermore, the spin-down rate before and after the giant flare is almost the same since the moment of inertia and the 
external magnetic fields are nearly constant. 

The first point that one may wonder is why the spin-down occurs with a release of energy. In a rotating star case, such as 
in the starquake model for the ordinary glitches (Ruderman 1969; Baym & Pines 1971), the energy release is always associated 
with the spin-up, since a rotating star has less oblate equilibrium shape as the star slows down. Nevertheless, it is not a trivial 
problem whether or not the spin-down occurs with a release of energy in a magnetized star case, since it is possible for a 
magnetized star to be prolate rather than oblate. Therefore, in this paper we shall address the following problem as a first 
step: can the spin-down occur with a release of energy in a magnetized star ? To approach this problem, we consider the most 
idealized model. We prepare several equilibrium polytropes with different magnetic configurations assuming axisymmetry and 
no rotations. Then, we investigate the relation of the energy and the moment of inertia between these equilibria to find at 
least one example in which the moment of inertia increases with a release of energy. 

In Section H, we will study the equilibrium of a magnetized polytropic star. In Section H, we will obtain the expression 
for the energy and the moment of inertia tensor. In Section W, we will compare equilibrium configurations and show that the 
spin-down with a release of energy is possible. We will also apply the results to the giant fiares of SGRs. Section pi is devoted 
to discussions and summary. We will consider observational signatures to test this model in the future. 



2 THE EQUILIBRIUM OF MAGNETIC POLYTROPES 

The equilibrium of a magnetized star has been studied by several authors since the pioneering work of Chandrasekhar and 
Fermi (1953). If one assumes a barytropic equation of state, one can make some progress in the understanding of the magnetic 
field on the equilibrium configurations. The equilibrium of a polytropic star with an axisymmetric poloidal magnetic field has 
been examined by Monaghan (1965, 1966a, b), and has been treated in a systematic manner by Trehan & Billings (1971). 
The structure of the magnetic field with a toroidal and a poloidal component in a polytropic star has been examined by 
Roxburgh (1966), and the equilibrium of a polytropic star with a toroidal and a poloidal magnetic field has been studied in a 
systematic manner by Trehan & Uberoi (1972). These models are constructed using a perturbation technique with the ratio 
of the magnetic energy to the gravitational energy as the perturbation parameter. This is essentially similar to that developed 
by Chandrasekhar (1933) and Chandrasekhar & Lebovitz (1962) for a slowly rotating polytrope, in which the perturbation 
parameter is the ratio of the rotational energy to the gravitational energy. Since the ratio & of the magnetic energy to the 
gravitational energy is small as long as H ^ W^^ {R/10^ cm)~''' (M / Mq)^ G, this method is the most suitable for our purposes. 
As we can see from equation (pi), it is not sufficient to consider only the first-order in a small parameter 5. Therefore we will 
extend and generalize the previous works to the second-order in a small parameter S. 

2.1 Basic equations 

The basic equations governing the hydrostatic equilibrium of a perfect conducting fiuid are 

= -Vp + pV<l>+-!-(Vxi3")xi3", (4) 

V^<E> = -47rGp, (5) 

V-if = 0, (6) 

where p is the pressure, p is the density, $ is the gravitational potential and H is the magnetic field. We assume the polytropic 
equations of state, 

' If the magnetic fields were well tangled, the deformation would not be elliptical but radial so that the energy shift is of order 5. In 
this case, our model gives too much energy. 
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p = Kp'+'^\ (7) 

and we normalize the quantities as 



p = pcO", r ^ a^ = 



K{n + l)p, 



-l + l/ri 



47rG 



1/2 



C, H = {A%GSf''^ pcah, $ = 47rGaVc0, (8) 



where 0, ^, /i and (j) are dimensionless quantities, n is the polytropic index, pc is the density at the center, and ii' is a constant. 
Then, equations (W), (Bl) and (q) are written as 

. _ve + v^ + 5i^2i^, (9) 

v^0 = -e", (10) 

Vh = 0, (11) 

where 5 represents, apart from a numerical factor, the ratio of magnetic to gravitational energy, and V is the divergence with 
respect to ^. 

We shall assume that the magnetic field is axisymmetric about the z-axis. Then we can express the magnetic field as 
a superposition of a poloidal and a toroidal field in terms of two scalar functions P{S,,p) and T{£,,p) in the manner (e.g. 
Chandrasekhar 1961), 

.^^^=Vx{mPe^)+mTe^, (12) 

where e^p is a unit vector along the (^-direction and {zij, (fi, z) denotes the cylindrical coordinate, which is related to the spherical 
coordinate {^,9,ip) as zu — ^sinS — ^(1 — p^Y''^ and z — ,^cos0 = ^p. Note that the projection of the lines of force on the 
meridional planes gives -cu^P = const. Since equation (H) requires Vx [{\7xh)xh/Q"] = 0, it can be shown that (Woltjer 
1959; Ferraro 1954; Parker 1979) 

= \/[e-(t>-SNp{zu'^P)], (13) 

AsP = -e"N'p{zu^P)-TNT{T^^P), (14) 

zu^T = Nt{zu^P), (15) 

where Nt{x) and Np{x) are arbitrary functions which characterize the configuration of the magnetic fields, and 

'^ ^ dzu'2 '^ vj dzo ^ dz-^ ^ dS,'^ ^ C 9C "*" C^ dp^ ^^ dp' ^ ' 

is the five-dimensional axisymmetric Laplacian. Taking the divergence of equation ( |l3| ) , we obtain 

\/^e^~e" + s\/'^Np{uj^p), (17) 

with equation (M). 

2.2 Perturbative approach 

We shall suppose that the magnetic field is so small that S may be treated as a small perturbation parameter. On this 
assumption, we will obtain the solutions for the density 0, the poloidal field P and the toroidal field T to the second-order 
in S. We assume that the solutions for O, P and T have the forms, 

e = Qo+Sei+ 5^02 + 0(5'-'), (18) 

P = Po + 5Pj + 0{5^), (19) 

T = To + STi + 0{5^). (20) 

By substituting these equations into the basic equations (|l4), (llq) and ([17|), we obtain in each order of S 

v'eo = -e;^, (21) 

V^Gi = -ner'ei + V^iVp(t;<7^Po), (22) 

v^ea = -ner'e2-^^^^^er^e? + v^[.^^Aiv;,(w^Po)], (23) 

A5P0 = -e^N'p{m^Po)~'ToN!r{ru''Po), (24) 

A5P1 = -nQo'^QiN'pizu'^Po) -Qozu^PiN'p{zu'^Po) -TiNTi^u'^Po) -Tozu^PiNt{vj^Po), (25) 

i:a2ro = Nt{uj^Po), (26) 

ca^Ti = t;<7Vi7Vr(i:i7^Po). (27) 
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These equations are to be solved order by order. It is obvious that depending upon the choice of the arbitrary functions Np (x) 
and Nt(x), there can be many solutions for the functions P and T which are consistent with the demand for equilibrium. 
Therefore it is convenient to make general discussions before specifying the functions Np(x) and Nt{x). So, we will consider 
the formal solutions for the density Q and the gravitational potential </> in the next section, and the formal expressions for 
the energy and the moment of inertia tensor of the magnetic polytrope in Section H, without solving the magnetic fields P 
and T. We will specify the magnetic fields P and T in Section U. 

2.3 The solution for B and (j) to second-order in S 

We will obtain the formal solutions for the density O and the gravitational potential (j) before considering the specific magnetic 
fields. Firstly we perform the Legendre expansion as 



oo 
oo 
oo 

zu^'PiN'pizu^Po) = ^r™(OP™(M), 



(28) 
(29) 
(30) 
(31) 



where Pm{lJ.) are Legendre polynomials of order m. Then, from equations (Ell), (0) and (E3), we find that the radial functions 
Qo, 'i/'m and 7m, satisfy the equations 



Vmi^pm - *m) = -nQo'^ipm, 

OO 

^m(7m— Fm) = — 720o~"^7m Bq'^ / ^j^lQh 



(32) 
(33) 

(34) 



where 

'" " e d^ V d^ 



^ + 1) 



e 



and Qim. is defined by 

oo oo 



(35) 



(36) 



Note that Qio = ^i/(2/ + 1) and Go in equation ( |32[ ) is, of course, the Lane-Emden function of index n. The radial functions 
00, 1pm and 7m in equations (p2|), (BSl) and (M) are to be solved by imposing the boundary conditions 



eo(o) = i, e;(o) = o, ^„(o) = o, ^Uo) = o, 7^(0) = o, 7^(0) = 0, 

since we shall let pc denote the central density. 

With equations (Eq), (p9|), (B(J) and (pl|), the solution for 4> can be obtained by integrating equation (113) 



= Co + eo + 5 



Cl;0 + ^{i>m- *m) PmilJ.) 



+ 5' 



C2;0 + ^ (7m - Tm) Pm{fl) 



(37) 



(38) 



where we write the integral constant as cpo = Co + 5Ci;o + 5^C2;0- The gravitational potential outside the polytrope satisfies 
Laplace's equation V^^ — 0, where the tilde in (f> distinguishes the exterior potential (j) from the interior one (ji. The solution 
of Laplace's equation which should be associated with the interior solution dSq) is, clearly 



Ko 



e 






= T^+^E^^-(/^)+^^E^^-(^)- 



K2.„ 



(39) 



It remains to determine the various constants which occur in the solutions for (p and by imposing on them the boundary 
conditions. 
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Ih(m) 



Ih(m) ' d^ 



h(m) 



ac 



(40) 



H(m) 



where the boundary of the polytrope H(/i) is given by 

oo oo 

H(m) = Co + (5 ^ q^Pmiti) +5^^ t„,P„,{fi). (41) 

m — Tn^O 

Here ^o is the first zero of the Lane-Emden function. The requirement that the density vanishes on the boundary surface 
0(H(/i),/i) — determines the constants qm and tm as 

V'm(Co) 



ei(Co) ' 



eo(Co) 



e;,(Co) 



Co 



Qim(?o) 



(42) 
(43) 



where Qim is defined in equation (JS^)- Then, evaluating <^ and (f> and their derivatives on H(/i), the boundary conditions in 
equation (|4C|) give 

^0 = -eoe;(eo), Co = -eoe;(Co), 
^i;™ = — ^C+'[^™(?o)-*™(Co)], 

m + 1 
Ci;o = -^o(Co) + *o(Co) - Co [V'o(Co) - *o(?o)] , 

if2;™ = ^C'[7™(?0)-r:„(C„) + Q„], 

C2;o = -7o(?o) + ro(Co) - Co [7o(eo) - r;(Co) + Qo] , 

and the additional boundary conditions for tprn and 7™ with m > 1, 

(m + 1) [^-^(50) - *^(Co)] + Co K(Co) - *:n(eo)] = for m > 1, 
(m + 1) [7„(Co) - r„(Co)] + Co [7,'n(eo) - r;„(Co) + Q,n] = for m > 1, 
where 



Q 



2e^(eo) 



■X]^'(?o)Q'-(?o), 



(44) 

(45) 
(46) 

(47) 



and Qim is defined in equation (pq). For rz < 1, this constant Qm diverges since 0o(Co) = 0, which means that the expansion 
with respect to S breaks down.rj Hereafter we will consider the case of n > 1. 

Now we get the formal solution of the problem, i.e., we can calculate the density Q and the gravitational potential <j} if 
the magnetic fields "^m and r,„ in equations (pOl) and (|3l|) are given. 



3 ENERGY AND MOMENT OF INERTIA TENSOR 

In this section, we will obtain the formal expression for the differences in the total energy and the moment of inertia tensor 
between two equilibrium configurations. The energy and the moment of inertia tensor will be obtained in the expanded forms 
with respect to the small parameter 5. 



3.1 Mass 

The total mass of the fluid is given by 

3/2 



M= pdr = 
'v 



K{n + 1) 



AtvG 



-l/2+3/2n 



e"d^e- 



(48) 



+ Chandrasekhar's perturbation technique breaks down near the surface of the polytrope since |©o| < |50i| at ^ ~ ^o- For a strict 
argument, some matching procedures may be necessary as in the study of a slowly rotating star (Monaghan & Roxburgh 1965; Smith 
1975; Singh & Singh 1984). However, from the study of a slowly rotating star, it seems that the perturbation technique gives a sufficiently 
accurate answer for, at least, integrated quantities such as mass (Anand 1968; James 1964). So we leave strict arguments in the future. 
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In the following discussions, we wish to compare configurations of the equal mass. Since the magnetic field causes the density 
to change, the equality in the masses can be achieved by an adjustment in the central density, 

pc = po + 5 pi + 5^P2- (49) 

With equations (tq), (^81), (M) and (ksl), the total mass in equation (Uq) can be expanded as 

M = Mo + 5Ah+5'^M2, (50) 

where 



Mn = 



Ml = 



Ma 



€o 



«o 



e^fd^, 




2n po 




m — 


^\i'rr.f 



2n po 2n 



P2 n + 3 f pi 



po -in \ Po 



The mass is normalized by the unit 

Cm = 47rpo«0) 
where 



ao 



47rG 



1/2 



(51) 
(52) 

A/o-?oQo. (53) 
(54) 
(55) 



is the unit of the length. Note that Qo in M2 comes from the upper limit in the integral, which is not ^0 but H(/^), i.e., the 
integral in equation (ksl) has to be performed as J 0"cfi£_ — 2tv J dp J" '^ d^^^S". The requirement of the equal mass is 
assured by Mi = and M2 = 0, which determine pi and p2 as 

2n i^'oi^o) - %{^o) 



Pi 



3-n e(,(eo) 



P2 _ 2n 7o(Co)-r;,(Co) + Qo _^n(n + 3) 



po 



3 -- n 



e;,(Co) 



(3-n)2 



i'oi^o) ~ %{^o) 



e^(Co) 



(56) 
(57) 



where we use equations (Al), (A3) and (A8) 



3.2 Energy 

The total energy of this system is given by (e.g., Chandrasekhar 1961; Woltjer 1959) 

S ^M+U + W, 

where 



M = 



StvCe 



u 



w 



Ce Jv 47r \ po 



2Ce J^'^ 47r I po 



47r \po J 

l/2+5/2n 



'V 

-1/2+5/2II 






(58) 

(59) 
(60) 
(61) 



are the magnetic energy, the internal energy, and the gravitational energy respectively, and we use equations (M), (pi) and ([13). 
The energy is measured in the unit of 

Ce = 47rA:(n + l)pj+'''"ag, (62) 

where ao is given by equation (pq). 



3.2.1 Magnetic energy 

With equations wM, (BOl) and (Ea), the total magnetic energy in equation (M) can be expanded as 

M =SMi + S^M2, 
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where 

M2 = ^ I (-^zu^PiAsPo - ^uj^PoAsPi + zu^TiTo) d'C + ^-^—Mi. 
47r y^ V 2 2 J 2n po 



By using equations (^6|), ( A21 ) and ( A37 ), the above expressions for A4i and A42 can be reduced to 

/•Co 



eo*o?'d^, 



M2 



«0 



3-n e(,(Co) 



J ° e'^^'^ed^ + / " E ^^^TTT®^"''^'"*'"^''^^- 



(64) 
(65) 

(66) 
(67) 



3.2.2 Internal energy 

With equations ([l8|), (|28|), ( |29|) and ([l9|), the internal energy in equation ( pCJ) can be expanded as 

U^Uo + 5Ui+ S'^W, 

where 



Wo = 



Wi 



n + 1 



io 



Q^'Cd^, 



«o 



2n po 



f „2 



rrfc + 



5 — n /9i 5 — n 

-Ml 



p2 n + 5 f p 



Po 4n y po 



Wo. 



2n po 2n 

By using equations {pm, (pTI), (|A2|), (A4) and (A9), the above expressions for Wo, Wi and W2 can be reduced to 

Wo = ^eo[e(,(Co)]', 



Wi = 



3 — 71 



Co 



W2 = - 



eo*o?'rfC, 



«0 



eo nrdc - 



«0 



n I /■•„„w,3,t 5-n Vo(Co) - *o(Co) / ^n,T,/^3_ 



eo*ordC 



+ 



3-nl /„ " "' "' 3-n 6o(Co) 

^_n ^ ^ 



(68) 

(69) 
(70) 
(71) 

(72) 
(73) 



(74) 



Therefore, from equations (pd), (pTj), (jTSl) and ([741), there are following relations between the internal energy and the magnetic 
energy, 



Wi = - 

W2 = - 
where 

M2 = : 



3 — n' 

n 

3 — 71 



Ml, 
7\/(2 +7V(2 



1 , /-So °° 

— ^0 K(Co) - *o(Co)] +j Y. 2(2m + l) ^°"'^'" t'^'" + ^^ ("^^ ~ *™)] ^''^^■ 



(75) 
(76) 

(77) 



Note that M2 does not contain the second-order quantities, such as jm and Fm in equations (ES) and (Bl|), but the non-linear 
terms of the first-order quantities, such as ipm and Vl'm in equations (128) and (pOJ). 



3.2.3 Gravitational energy 

With equations (tlq), (^81), (p9|), (pq) and (M), the gravitational energy in equation ( pi] ) can be expanded as 

W = Wo + (5Wi + 5^W2, (78) 
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where 



Wo 
Wi 

W2 



Co 



Co 



nej}- Vo (Co + 60) + eS (Ci;o + Vo - *o)] c'rfC + ^^— Wo, 



2n po 



(79) 
(80) 



Co 



nOo So + ^ 



"("-^)-er(^„)^ 



2(2m+l) 



I i'^ n — 1 



00 

^0Cl;0 + > —ri>m (Tpm ^ '^ r. 

^-^ 2m + 1 



+ izZ^^Wi+^-" 



2n Po 



2n 



P2 n + 5 f pi 



po 4n y Po 



(Co + 60) 

+ ef'?(C2;o + 7o-ro)U'rfC 



Wo + ico^gQo. 



(81) 



Note that Qo in W2 co mes f rom t he upper hmit in the integral, which is not ^0 but H(p). By using equations (|4j), ([Sq), (p7|), 
(A1)-(A5), (A^)-(AIO) and (A13), the above expressions for Wo, Wi and W2 can be reduced to 



S n 

Wo = Wo, Wi = Ml, W2 

n 3 — n 



X2 + 7; M2- 



3 — n 



3 — n 



(82) 



5.^.-</ Total energy 

For configurations in static equilibria, the virial theorem for the hydromagnetics is given by (Chandrasekhar 1961) 



W + M + -U ^ 0. 

n 



(83) 



We have shown this theorem explicitly to the second-order in 5 with equations (iTq) , (iTq) and (K2) . By eliminating W + M in 
equation (|58), the total energy can be expressed as 



3 — n 



U 



- — -Uo + SMi+5^ {M2+M2 



(84) 



with equations ([75|) and WW. This equation (pj) can be interpreted as that the total energy is the sum of the zeroth-order 
background energy £0 = —[(3 — n)/n]Uo and the total magnetic energy M = 5Mi + S^M2 with a correction S^M2 which 
comes from the product of the first-order quantities. 



3.3 Moment of inertia tensor 

The moment of inertia tensor is diagonal in the chosen representation. In view of the axisymmetry, the two distinct components 
of this tensor areP 



Xii 



1 / 1 / \ -3/2+5/2" 

T ^ / 2/, , 2n ,3 J^ ( Pc \ / ^n,2/, , 2n ,3/ 

I22 = 7:7;- I pr [1 + fi )d r - ' I Q ^ /I , .. \j ^ 



2C/ 



stt y po 



e"e'(i + /^')d'C, 



, /- 1 / \ -3/2+5/2n ^ 

X33 ^ i//^^(i-M^)'^^-i(S) Xe"e^(i-^Vc, 

where the moment of inertia tensor is measured in the unit of 

C'l — 4npoao. 

With equations (tlq), (tsh, (h^) and (ka), we can expand the moment of inertia tensor as 

111 ~ T^Q + <5Xll;l, 

2^33 ~ I0 + 5X2,3-1, 



(85) 
(86) 

(87) 



(89) 



8 The moment of inertia tensor can be constructed from the second moment of the mass distribution lij by Xij = Sijtr(Iab) — hj- The 
i-component of the angular momentum is given by Ji = lij^j where f2j is the j-component of the angular velocity of the rotation. (Of 
course, our analysis neglects the rotation.) 
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where 

9 r^° 



Ill;l 



I33;: 



«0 



Co 






5 - 3n pi 

5 - 3n pi 
2n Po 



(90) 
(91) 
(92) 



As we will see later, it is sufficient for our discussions to expand the moment of inertia tensor to the first-order in S. 



3.4 Differences in energy and moment of inertia tensor between equilibria 

We shall obtain the formal expressions for the differences in energy and moment of inertia tensor between two equilibria which 
have different configurations of the magnetic fields. In order to compare two equilibria, it is not sufficient to fix the mass. 
We need two relations of parameters between equilibria because there are two free parameters, the central density po and the 
ratio of the magnetic energy to the gravitational energy S, in the magnetized polytropic model. We shall take the magnetic 
energy as one more fixed parameter. Let us characterize the configuration of the magnetic field by A, and consider that the 
magnetic configuration changes from Xi to A/ . Hereafter let the subscript i (/) mean the initial (final) state. If we assume that 
the magnetic energy decays by some amount AA4 in this transition, the relation of the magnetic energy between equilibria is 
given by 



M{Xf) -M(\^) = [SfMl(\f) + SJM2{\f)] - [S^Ml{\^) + S^ M2iX^)] 



-AM, 



(93) 



which determines the relation between Si and 5/ . Since the amount of the decaying magnetic energy depends on the process 
of the magnetic reconfiguration, we will treat AM as a parameter in our discussions. Then we find from equations (184) and 
(BSn that the difference in the total energy between two configurations A; and Xf is given by 



A£_ ^ £{Xf)-£{X^) ^ iMliM 
\£o\ ~ \£o\ ' ' \£o\ 



M2{Xf) M2{X,) 



AM 



MliXf) MUX,) S^MliX,) 



(94) 



where £0 is the lowest order of the total energy, £0 = —[(3 — n)/n]Uo, and we keep the lowest order of terms which include 
AA^. The energy is released when A£ < 0. If we can neglect the term of AA^, i.e., | AA1| ^ S , the released energy is of order 
lAf/fol ~ 5^. 

From equations (pq), ( p9| ) and d93), the difference in the moment of inertia tensor between two configurations Xi and A/ 
is given by 



AXkk _ IkkjXf) —XkkjXi) _ ^ MljXi) 



1hk;l{Xf) Jfefc;l(Ai 



AM 



l(A/) 



Xi(A/) Mi(AO 5,Mi{X,) Mi{Xf) 



(95) 



keeping the lowest order in 5. The moment of inertia tensor increases when AXkk > 0. If we can neglect the term of AM, i.e., 
jAA^I ^ 5, the difference in the moment of inertia tensor between two equilibria is of order \AXkk/XQ\ ~ 5. 

From equations (M) and (p)5[), the conditions for the moment of inertia tensor to increase with a release of energy, A£ < 
and AXkk > 0, are reduced to 



M2{K 
MUX, 



+ 



V 



S,Mi{X,) 

1 2'fcfc;l(Ai) 

^ Xi(AO 



> 



< 



M2{Xf) 
MliXfY 

Xkk;l(Xf) 

MiiXf) ' 



where 



V 



AM 



<1, 



(96) 
(97) 

(98) 



5iMi{Xi 

is the ratio of the decaying magnetic energy to the initial total magnetic energy. Note that these conditions depend on the 
zeroth- and the first-order quantities and not on the second-order quantities, i.e.. Mi, M2 and Xkk-.i do not depend on 7m 
and Fm but on &o, ij-'m and ^m as we can see from equations (pq), (|77]), (|9l| ) and (P2[). 



4 COMPARISON OF EQUILIBRIUM CONFIGURATIONS 

So far we have obtained the general solutions and expressions for the equilibrium configurations without solving magnetic 
fields. In this section, we specify magnetic configurations to investigate the conditions for the moment of inertia tensor to 
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increase with a release of energy in equations ^61) and (pTf). We shall obtain the magnetic fields to the lowest order in S, i.e., 
Po and To, since the conditions in equation (Bq) and (b7[) do not depend on the second-order quantities as noted in the last 
of the previous section. 



4.1 The solution for P and T to first-order in 5 

We shall consider a special case of the magnetic fields characterized by 

Np{x) = -X, Nrix) = Aa;, (99) 

in which case the solutions for Po and To are functions of only the radial coordinate ^ and are easily calculated (Woltjer 
1960; Trehan & Uberoi 1972). Then A controls the configurations of the magnetic fields. Substituting the above expressions 
for Np{x) and Nt{x) to equations (M) and (bq), we get equations to be solved, 

AsPo + A'Po^e^, To = APo, (100) 

where Oo is the Lane-Emden function of index n. Since we are considering the larger interior fields compared with the exterior 
fields, the most convenient way to impose the boundary condition is to require that the fields vanish at the boundary (Trehan 
& Uberoi 1972), that is to the lowest order in 5, 

Poi^o)=0, ^(Co) = 0. (101) 

Then, the solution for Po is given by (Woltjer 1960; Trehan & Uberoi 1972) 

MO = tM^O eo{x)jiiXx)x^dx+-ji{X^) ej;{x)ni{Xx)x^dx, (102) 

where A is a root of 

eoix)ji{Xx)x^dx^O, (103) 





and jm (0 and Um (^) are spherical Bessel and Neumann functions of order n respectively. The first few roots of equation ( lOJ ) 
are given in Table 1 for the poly tropic indices n = 1, 1.5, 2 and 2.5. 

Having determined the solution for Po, we can calculate ^m in equation (|3C|) by 

*o = -|?'Po, *2 = |c'^o, *m = for m/0,2. (104) 



Then the radial functions ipm can be calculated by solving equation (p3) subject to the boundary conditions in equations (37) 
and (15). The quantities, Mi, M2 and Ikk-.i, can be calculated|j from Oo, 4'm and 'I'm by using equations (p6[), ([7 7]), (91) 
and (^. 

Now, we are in a position to examine the conditions for the moment of inertia tensor to increase with a release of energy 
in equations (^6|) and (|97|). 



4.2 Results 

For the moment, we neglect the decaying magnetic energy, i.e., r; = in equation (|9q). Note that ?; is a free parameter in our 
analysis since we need two relations of parameters (here the mass and the magnetic energy) between equilibria. We evaluate 
AI2/AI1 (energy) and Tkk;i/Mi (moment of inertia) by taking various values of A to examine the conditions for the moment 
of inertia tensor to increase with a release of energy in equations ( |96[ ) and (|97|). In Fig. 1, we plot Al2/A^i and Ikka/ Mi as a 
function of A for the cases of n = 1 and n = 2.5. A takes discrete values in Table 1. We can see that the conditions in equations 
( Pq ) and (BTI) are satisfied for the moment of inertia tensor of 33-component, I-yj,-!, i.e., if the magnetic configuration changes 
from the state of larger A to that of smaller A, M2/ Mi (energy) decreases while l-iz-^i/Mi (moment of inertia) increases. This 
means that the total energy decreases while the moment of inertia tensor along the magnetic axis increases from equations 



(M) and (|95|) . Therefore the answer for the question posed in Section 1.2 is that the spin-down with a release of energy is 
possible in a magnetized star. 

™ In the numerical calculation, we can solve all differential equations simultaneously using the fifth-order Runge-Kutta method (Press 
et al. 1992). 
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Note that the configuration of the polytrope is prolate rather than oblate since the moment of inertia tensor of 11- 
component Xii;i is larger than that of 33-component Xss;! from Fig. 1. Therefore, as the total energy decreases, the configu- 
ration of the polytrope becomes more spherical. This may meet our intuitions from the view that the self-gravitating fluid is 
the most stable at the spherical configuration. 

In Fig. 2, the projection of the lines of force on the meridional plane, which is given by S^'^zu'^Po ~ const, is shown for the 
case of n = 1 (see equations (4.17) and (6.17) of Parker 1979). The lines of force lie in the two dimensional surfaces which are 
made by rotating each loop in Fig. 2 about the symmetrical axis. For A = A2 the inner projected lines rotate in the different 
direction to the outer ones on the meridional plane. The lines of force for larger A have the more complicated structure than 
that for smaller A, i.e., the number of the smaller structure like the inner loops for A = A2 increases as A increases. Note that 
the flux of lines of force through any slice of the star is always zero because of the boundary condition /i(H(/^), ^) — 0. 

To evaluate physical quantities, we adopt the canonical values, Cm Mo = IAMq and R — Qo^o = 10^ cm, for the total 
mass and the radius of the polytrope respectively, which determine po and K from equations (m) and (Ba). Furthermore we 
adopt AT33/T0 = 10~* since this value is required for the positive period increase of SGR 1900-1-14. Then we can calculate 
Si from equation (Ba), and hence the released energy Ce^£ from equation (b4), where Ce is the unit of the energy in 
equation (p2). In Fig. 3, we plot the released energy CeAf as a function of the various sets of the initial and final state, 
(Ai, A/) = 10(/ - 1) - /(/ - l)/2 + i- f, where / = 1, • • ■ , 10, i = / + 1, ■ • ■ , 10 and / < i (see Table 2). We can also 
obtain the magnetic energy M — SMi required for AT33/X0 = 10~* from equation (psf). In Fig. 4, we plot the typical 
magnetic field H, which is defined by (//^/87r)(47ri?^/3) = CeM, as a function of the various sets of the initial and final 
state, (Ai, A/) = 10(/ — 1) — /(/ — l)/2 + i — f (see Table 2). These results in Fig. 3 and 4 are essentially the same as the 
order-of-magnitude argument in Section |l.2| . Note that there are relations, CbAiS oc (AXss/Xq)^ and H ex (AX33/Xo)^ • 

When the decaying magnetic energy is not zero, r; 7^ 0, it can occur that the conditions for the moment of inertia 
tensor to increase with a release of energy in equations (661) and (Wn) are not satisfied. In Fig. 5, we obtain the region 
of r; where the conditions in equations (bq) and ( |97] ) are satisfied as a function of the various sets of the initial and final 
state, (Ai, A/) = 10(/ - 1) - /(/ - l)/2 + i- f, where / = 1, ■ ■ • , 10, i = / -Hi, • • • , 10 and / < i (see Table 2). Here 
we assume AX33/ZD ~ 10^^ and consider the 33-component of the moment of inertia tensor. From this Figure, we find 
that the conditions in equations (M) and (pq) are satisfied when —10"^ ^ r; ^ 1 for n — 1 and —10^^ ^ ?? ^ 10^"^ for 
n = 2.5. Note that the decaying magnetic energy is comparable with the released energy when \ri\ ~ 10~* ~ AT33/T0 since 
|AA4| ~ \V-^\ ~ ^\{-M/^)£\ ~ "^^l^l- Therefore the moment of inertia increases with a release of energy if the decaying 
magnetic energy is less than the released energy. 

The numerical values of the various parameters are given in Tables 3 and 4. We can calculate such as A£ and Alkk/Io 
using these Tables and equations (M) and (|9q). We also give the ratio of the energies of the toroidal and the poloidal 
components of the magnetic field, A4Ti/Mipi. To the lowest order in S, it is given by 

1 / /I 5„2\ .3 i A 



4^ „. x^ . . ,0 



(i^^To^) d'C = Y / Pk'd^, (105) 



Mpi = ^ [ (-it^^'PoAsPo) d^'C = Mti -\ r QoPo^dC (106) 

From Table 4, the polytrope becomes more spherical as the ratio A4ti/Mpi decreases, which is a consequence of the fact 
that toroidal fields tend to make the star prolate and poloidal fields tend to make it oblate. It is interesting to note that the 
internal magnetic field generated by a post-collapse a-Q dynamo is probably dominated by a toroidal component (Duncan & 
Thompson 1992; Thompson & Duncan 1993). 



5 DISCUSSIONS 

In this paper, we have proposed one possible mechanism for the giant flares of the SGRs within the framework of magnetar, 
motivated by the positive period increase associated with the August 27 event. We assume that the global reconfiguration of 
the internal magnetic field suddenly occurs. Then, the shape of the magnetar is deformed and this deformation causes the 
moment of inertia and the total energy to change. With the internal magnetic field oi H ^ lO^'^ G, we can explain the positive 
period increase APt/Pt ~ 10"** as well as the energy ^ 10'*'* ergs of the giant flare of SGR 1900-1-14 if the magnetic fleld does 
not decay so much in this giant flare. In this mechanism, the energy source is not the magnetic energy but the gravitational 
energy. In order to investigate whether or not the spin-down occurs with a release of the energy in a magnetized star, we have 
analyzed the magnetic polytrope as the most idealized model using the second-order perturbation technique. We have found 
that it is possible for the moment of inertia along the axisymmetric axis to increase with a release of energy. This result can be 
understood as follows. The polytrope is prolate at the first state. As the polytrope becomes more spherical, the total energy 
decreases since the spherical configuration is the most stable. On the other hand, the moment of inertia increases since the 
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radius in the equatorial plane increases. Hence the spin-down is a natural consequence of the internal field reconfiguration. It 
is interesting to note a suggestion about prolate white dwarfs (Katz f989). 

In order to understand what is physically going on in our model, we would like to describe typical physical quantities. 
Let us consider a star with a polytropic index n = 1, a mass Cm Mo — IAMq, and a radius R — ao^o = 10^ cm. The central 
density is given by po ~ 2.f9 x fO^^'g cm" from equations (M) and (pq), and the mass conservation is always assured by 
equations ( p6| ) and (pTf). We shall consider the first and final state of the magnetic field which is characterized by A = A2 and 
A = Ai respectively. The lines of force are shown in Fig. 2, and the fields vanish at the surface of the star. For simplicity we 
assume that the decaying magnetic energy is much smaller than the released energy, i.e., rj — AM = in equations (BSl) and 
(Pq). We shall adopt that the increase in the moment of inertia is ATss/Tq = 10~* since this value is required for the positive 
period increase of SGR 1900+14. Then, in the first state. Si = 1.89 x 10"'^ from equation (ba) and Tables 1 and 4, the central 
density is po + 5pi — po + 1.58 x lO^^g cm""^ from equations (|49) and (pq), and the star shape is prolate since the moment of 
inertia tensor of 11-component is larger than that of 33-component, T11/T33 = 1 + 3.25 x 10~*, from equations (BSl) and ( pEJ ) and 
Tables 1 and 4. In the final state, 5/ = 4.46 x 10~*, the central density is po + 5p\ = po — 2.05 x lO'^^g cm~'^, and the star shape 
is prolate since T11/T33 = 1 + 1.31 x 10~*. In this magnetic reconfiguration, the star becomes more spherical, the magnetic 
energy CeM — 9.60 x lO^^ergs is assumed to be nearly conserved, and the released energy is Ce\A£\ — 1.73 x lO^^ergs from 
equation (|94) and Table 4. 

The greatest difficulty in our model is the assumption of the sudden internal field reconfiguration. Thompson & Duncan 
(1995) considered that the internal field may diffuse into a configuration where an interchange instability is no longer inhibited 
by topology (Flowers & Ruderman 1977), since the diffusion of the strong magnetic field will be efficient (Goldreich & 
Reisenegger 1992; Thompson & Duncan 1996). We may speculate other possibilities, such as magnetic reconnection (Thompson 
& Duncan 1995, 1996), unpinning of the magnetic vortices in analogy to the terrestrial hard superconductors (e.g., Sauls 1989), 
and so on. In any cases, there are at least two theoretical problems. First, we have found that the decaying magnetic energy 
during the reconfiguration has to be less than the burst energy AM ^ 10*^ ergs, i.e., the magnetic energy has to be nearly 



conserved |AA^/A^| ;$ 10~ (see Section 4.2). Second, it is also necessary for the frozen-in condition to break between 
equilibria since each equilibrium has different Np{x) and Nt{x) from our analysis (Woltjer 1959). These constraints may be 
too strong to give up our model, or our model may be too simple. However, the interior of the neutron star is not yet well 
understood so that at present we can not say definitely whether such magnetic reconfiguration can occur or not. Therefore 
we had better await the future observational tests (see below) before reaching the decision. Anyways, it is an important fact 
that there exist different neutron star states separated by the right energy and the right shift in the moment of inertia. 

The first observational signature of our model is a pulsation of ~ ms period in the burst profile of the giant fiare. In our 
model, the released energy will be firstly converted into the oscillational energy of the neutron star (Ramaty et al. 1980). 
Since the damping time of the oscillation is about ~ 0.1 s (see below), a pulsation of ~ ms period in the initial pulse of the 
burst profile is expected. Indeed, there is the pulsating component with ~ 23 ms period in the hard initial pulse of the March 
5 event, which may be the observational signature of the oscillating neutron star (Barat et al. 1983; Duncan 1998). In the 
August 27 event, the data on the initial pulse with the fine time resolution are not available (Hurley et al. 1999; Feroci et al. 
1999; Mazets et al. 1999). Since the oscillations may give us the evidence for the superstrong magnetic field (Duncan 1998), 
it is important to observe the giant flares, especially the initial pulse, with the fine time resolution. 

The second observational signature of our model is the gravitational wave. According to our model, the released energy 
is related to the magnitude of the period increase. The predicted energy <: 10^^ ergs in Fig. 3 is somewhat larger than 
the observed energy ^ 10*^ ergs of the gamma-ray. However, it is interesting to note that the observation of the Earth's 
ionosphere indicates the presence of an intense initial low energy photon which carries substantially higher (by a factor of 
~ 9) total energy. Furthermore, the nonradial (p- and /-mode) oscillations, which will be excited in our model, will be 
principally damped by the gravitational waves|ll| (Ramaty 1980; McDermott, Van Horn & Hansen 1988; Cutler, Lindblom & 
Splinter 1990). Decay of these oscillations by gravitational waves may account for the duration of the initial pulse, ~ 0.1 
s (Ramaty 1980; Lindblom & Detweiler 1983). The characteristic gravitational wave amplitude (Thorne 1987) is given by 
he = h^/n — h^fjr = {4,GFt/-kc? fY''^ where / is the gravitational frequency, t is the damping time of the oscillation, n is 
the number of observed cycles and F is the gravitational wave fiux. Since the gravitational wave fiux can be estimated from 
the available energy E and the distance d of the source \yy F — E/A-Kd? — E /^nd^T, we obtain 

I We can make a rather crude estimate as follows. The luminosity of the gravitational waves can be estimated by the quadrupole formula. 
Low = (G/5c5) /fjkfjk) ~ {G/c'^) (^Af/AT^Y ~ 5^ (c^/g) {GM/c^FI^^, where Af ^SMB? is the quadrupole of the neutron star. 



and AT ~ yj B? jGM is the free fall time. We estimate the luminosity of the electromagnetic waves by the dipole formula. Since there 
are relations, I^H ~ const, p ~ HR^, and AR ^ 5R, where p is the magnetic dipole moment, we can obtain Lem = (2/3c'') |/i| ~ 

(l/c^) (Ap/AT^) ~ 5'^(cH'^R^) (GM/c'^R) . The ratio Lem/Lgw ~ S{c^R/GM) is small for typical parameters. Since the above 
estimate is too crude, we need a more precise calculation for the energy conversion efficiency from the neutron star oscillation to the 
gamma-ray. 

© 2001 RAS, MNRAS 000, 0-|^ 



14 K. loka 

'^-'^ ''"'" ( lo^^ii) ^'' (24) "'' (5]^) " • (^0^) 

Therefore, if the energy that is released as the gravitational waves is larger than ~ 10''^ ergs, which may be possible in our 
model from Fig. 3, the gravitational waves from SGRs within ~ 10 kpc could be detected (de Freitas Pacheco 1997; Mosquera 
Cuesta et al. 1998) by the planned interferometers such as LIGO (Abramovici et al. 1992), VIRGO (Bradaschia et al. 1990) 
and LCGT (Kuroda 2000). The event rate will be about 2/20 ~ 10"^events/yr. 

Although the energy source of the giant flares is the gravitational energy in our model, this energy may be supplied by 
the decaying magnetic field before the giant flares (Katz 1982). The occurrence of two giant flares in different sources within 
twenty years most likely suggests that such events occur in most SGRs once every ~ 100 years (Hurley 1999). Therefore, if 
the SGRs are active for ~ 10^ years, the total released energy as the giant flares is about ~ 10^*(-H'/10^®G)* ergs. This energy 
can be supplied by the total magnetic energy ~ 10*^(-ff/10^^G)'^ ergs. 

We have implicitly assumed that the magnetic axis is parallel to the spin axis. If the magnetic axis is inclined to the spin 
axis, the spin-up (not spin-down) may also occur at the giant flares since the moment of inertia along one axis within the 
equator decreases as the energy decreases from Fig. 1. However, there may be some tendency toward an initial alignment of 
the magnetic and spin axis if the magnetic field is generated by a post-collapse a-Q dynamcrj (Thompson & Duncan 1993; 
Thompson et al. 2000) . The electromagnetic (also gravitational) torque which brakes the rotation of a neutron star also tends 
to align the magnetic axis with the rotational one on the time scale for rotational braking (Goldreich 1970; Cutler & Jones 
2001; Melatos 2000). Furthermore, the angular velocity Q decreases (i.e., spin-down) as long as the inclination angle is less 
than 46 degrees when the polytropic index is n = 1 and the state changes from A2 to Ai , for example, since the change of T33 
is comparable to that of In from Table 4. In either case, the magnitude of the period jump is related with the released energy 
through equations (H) and (pi). The misalignment of the spin and the magnetic axis may also give rise to the modulation of 
the spin-down history due to free precession and radiative precession (Melatos 1999; Thompson et al. 2000), which would 
provide a direct measure of the internal magnetic field. 
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Al Integrals of So 

With equation (p3) and the boundary condition Go(5o) ~ 0, we can show that (e.g., Chandrasekhar & Lebovitz 1962) 



Co 



e^+^e^dc = ^Cone^(eo)] 



5 — n 



(Al) 
(A2) 



A2 Integrals of ^o and 'i'o 

We will show the following relations, 



io 



3 — 71 



7v, = / eSVor^C = T^ <^ dQ'oi^o) [^o(Co) - *o(Co)] - / QoKfd^ \ , 



?0 



Ho 



e[?*oC'rfC = CoeaV>o(eo)-*o(eo)] + (l-n)7^. 



(A3) 
(A4) 
(A5) 



Firstly, we can verify equation (A3) by making use of equation (|3^) satisfied by ipo. Secondly, we can verify equation ( A4) by 
performing partial integrals several times as 



3 , "(eo"Vo)'e'rfC = ^ 



" ner'Q'o^^oed^ - / " Qoi'ofd^ 



Co 



= ^Uoe[,(Co) [V'o(Co)-*o(Co)] - T" (ceo)'(V'o-*o)?'dc+ T" (C'e[,)'VoerfO 

= ll eoe[,(Co) [v^o(Co) - *o(Co)] + r Ooii^'o- *o) C'rf? + / " {2e'o + ^o'o) *oe'de 

= 1 1 eoe;(5o) [Vo(Co) - *o(Co)] - y^° Oo [e {i^'o - *(,)] ' d^ + f" {eOoY *oCdC i 

= l\ eoe;,(Co) [V'o(Co) - *o(Co)] +ni^~ r e^^fdd, 



(A6) 



with equations (p3) and (63), and the boundary condition Go(^o) = 0. Finally, we can verify equation (A5) by performing the 
above integration in a different way, 

nl^ = - / " Oo [e (V'o - *;)] 'd^= ["q'o {i'o - *;) fd^ = Coe^ [V'o(Co) - *o(Co)] + h ^ f " e^^o^'d^ (A7) 

Jo JO Jo 



A3 Integrals of 70 and Fq 

We will show the following relations. 



/■Co ria °° 

/ nerSoC'rfC = -eo[7o(Co)-r^(Co)]- / ^ 

"'0 Jo ^_n 
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rCo 



(0 
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I Y. 2(2m+l) ^°'''^'" t(^ ^ ''^'^'" + ^^'^™] ^'''^ 
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Qo^otd^ = eoO;, [7o(?o) - Fo(Co)] + (1 - n)/^ - 

JO 



«0 



E 



"("-^\er^(^^)^C^d5 



2(2m + l) 



(A8) 



(A9) 
(AlO) 
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Firstly, we can verify equation ( |A8[ ) by making use of equation (M) satisfied by 70 witli Qio = ipi/{2l + 1). Secondly, we can 
verify equation (A£) by performing partial integrals several times as 



^ 3 



€0 
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"'0 Jo Jo „_n 
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(All) 



with equations (p2|), ( |34[ ) and (|47|), Qio = ^i/i'il + 1) and the boundary condition 0o(^o) = 0. Finally, we can verify equation 
(AlO) by performing the above integration in a different way, 
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A4 Integrals of the product of tpm and ^,i 

We will show the following relation. 
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which is used to show equation (fed). By using equation (p3), the above equation (A13) is shown to be equivalent to 
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= ko K(eo) - *o(Co)] {M^o) - *o(Co) + Co K(Co) - *o(eo)] } 
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By performing partial integrals, we can show that 
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Using the above relations with the boundary conditions for tprn in equation (|45[), we find that equation (A14) is satisfied, and 
hence equation (A13) is satisfied. 



A5 Integrals of Pq and To 

We will show the following relation. 

First, the left hand side of the above equation can be transformed as 
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where we use equations 
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j) and (|30|). Noting that 
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and C(9/9C)a' — ■co[d/dzo)z + z^d/dz),^, we divide equation ( A22 ) into six parts, 
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(A22) 
(A23) 
(A24) 

(A25) 
(A26) 
(A27) 
(A28) 
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/mi-5 = - zu'^To^^ (zu^To) dzudzd(p, (A29) 
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By performing the partial integrals, we can show that 
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Then, the sum, {4-k) ^ X/m=o ^Afi;m, gives the right hand side of equation ( A21 ), and hence we have shown equation ( A21 ) 



A6 Integrals of Pi and Ti 

We will show the following relation 
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First, the first term in the left hand side of the above equation can be transformed as 
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where we use equations (p4)-(E7|) and (Blf), and 
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In a similar manner to the previous section |A5| , by performing the partial integrals, we can show that 
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zu TiTod ^. 



With equations (ESl) and (BQ), we perform the angular integrations in equation ( A41 ) to find 
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By combining equations (A38), (A.42), ( A43) and (A44), the proof of equation (A37) is completed 
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Figure 1. The quantities, M2/ M\, Xii;i/yVli and Izs-^i/ Mi, in equations ( |96| ) and ( p7| ) are plotted as a function of A for the polytropic 
indices n = 1 and n = 2.5. A controls the configuration of the magnetic fields and takes the discrete values in Table 1. As A decreases, 
A4.2/-M\ decreases, which means that the energy decreases, In^i/yVfi decreases, which means that the moment of inertia of 11-component 
decreases, and IzZ;! / -Ml increases, which means that the moment of inertia of 33-componcnt increases. 





Figure 2. The projection of the lines of force on the meridional plane, which is given by 5^' vj Pf^ = const, is shown for the case of 
n = 1. A controls the configuration of the magnetic fields and takes the discrete values in Table 1. As we can see from this Figure, the 
lines of force for larger A have the more complicated structure than that for smaller A. 

This paper has been produced using the Royal Astronomical Society/Blackwell Science lOTfTJX style file. 
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where / = 1, ■ • ■ , 10, i = / + 1, • • • , 10 and f < i (see Table 2), is shown. Here we assume AX33/X0 = 10"'*, which is required for the 
positive period increase of SGR 1900+14. 
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Figure 4. The typical magnetic field H, which is defined by (H^ /8it)(4:1tR^ /3) = CeM, as a function of the various sets of the initial 
and final state, (A;, \f) = 10(/ - 1) - /(/ - l)/2 + i- f, where / = 1, ■ • ■ , 10, i = / + 1, • • ■ , 10 and / < i (see Table 2), is shown. Here 
we assume AX33/X0 = 10~^, which is required for the positive period increase of SGR 1900+14. 
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5.510036 


A2 


3.407865 


3.183636 


4.298552 


5.988518 


A3 


4.430077 


4.053194 


5.047380 


6.611395 


A4 


5.443462 


4.915842 


5.765941 


7.187845 


As 


6.452475 


5.778637 


6.488668 


7.776904 


As 


7.458980 


6.641027 


7.210835 


8.362584 


A7 


8.463904 


7.503016 


7.933020 


8.949183 


Ag 


9.467764 


8.364664 


8.655152 


9.535587 


Ag 


10.47087 


9.226033 


9.377241 


10.12207 


Aio 


11.47343 


10.08718 


10.09929 


10.70856 



Table 1. The first few roots of A for the polytropic indices n = 1, 1.5, 2 and 2.5. 



(Ai,A/) 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


i 


2 


3 


4 


5 


6 


7 


8 


9 


10 


3 


4 


5 


6 


7 


8 


9 


10 


4 


5 


6 


7 


8 


9 


f 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


2 


2 


2 


2 


2 


2 


2 


3 


3 


3 


3 


3 


3 


(A«,A^) 


24 


25 


26 


27 


28 


29 


30 


31 


32 


33 


34 


35 


36 


37 


38 


39 


40 


41 


42 


43 


44 


45 




i 


10 


5 


6 


7 


8 


9 


10 


6 


7 


8 


9 


10 


7 


8 


9 


10 


8 


9 


10 


9 


10 


10 




f 


3 


4 


4 


4 


4 


4 


4 


5 


5 


5 


5 


5 


6 


6 


6 


6 


7 


7 


7 


8 


8 


9 





Table 2. The labeUing scheme of the x-axis in Fig. 3, 4 and 5, i.e., (A,;, A/) = 10(/ - 1) - /(/ - l)/2 + « - / where / = 1, • ■ • , 10, i = 
/ + 1, ■ • ■ , 10 and / < i, is explicitly shown. A controls the configuration of the magnetic fields and takes the discrete values in Table 1. 
The subscript i (/) means the initial (final) state. 
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24 K. loka 



1.5 



: 2.5 



«0 


3.14159 


3.65375 


4.35287 


5.35528 


Mo 


3.14159 


2.71406 


2.41105 


2.18720 


1^0 1 


1.57080 


8.64015 - 1 


4.45157- 1 


1.78659 - 1 


lo 


8.10448 


7.41315 


7.07403 


7.01315 



Table 3. A tabic of constants for the polytropic indices n = 1, 1.5, 2 and 2.5. §o is the first zero of the Lane-Emden function. Mq is 
the normalized mass of the zeroth-order. I^ol is the normalized total energy of the zeroth-order. Xq is the normalized moment of inertia 
tensor of the zeroth-order. The number following the plus or the minus sign indicates the power of 10 by which the table entry should 
be multiplied. This notation is followed in all tables. 







Ml 




M2/MI 


Xiiii/Xi 


2^33;l/A^l 


Mti/Mpi 


n = 1 


Ai 


1.30707 




1.57386 - 1 


4.09418 


2.27235 


1.86290 




A2 


3.07662 - 


- 1 


4.66603 - 1 


5.39981 


8.81646 - 1 


2.94224 




A3 


1.32171 - 


- 1 


6.55343 - 1 


6.18915 


-3.14410-1 


4.37118 




A4 


7.34761 - 


■2 


7.73213 - 1 


6.70086 


-1.18614 


6.15452 




As 


4.70193 - 


■2 


8.49530 - 1 


7.04562 


-1.80598 


8.29340 




Ae 


3.28329 - 


■2 


9.00955 - 1 


7.28585 


-2.25076 


1.07882 + 1 




At 


2.43167 - 


■2 


9.36929 - 1 


7.45842 


-2.57605 


1.36392 + 1 




As 


1.87852 - 


-2 


9.62941 - 1 


7.58583 


-2.81905 


1.68463 + 1 




Ag 


1.49784 - 


-2 


9.82294 - 1 


7.68219 


-3.00433 


2.04097 + 1 




Aio 


1.22402 - 


■2 


9.97052 - 1 


7.75663 


-3.14830 


2.43293 + 1 


n = 1.5 


Ai 


5.64307 - 


- 1 


4.77567 - 1 


1.02428 + 1 


6.25838 


2.42892 




A2 


1.77898 - 


' 1 


1.00202 


1.23688 + 1 


4.49438 


4.66722 




A3 


8.48260 - 


-2 


1.17620 


1.32906 + 1 


3.15561 


7.78376 




A4 


5.08044 - 


-2 


1.31401 


1.37993 + 1 


2.36568 


1.18083 + 1 




As 


3.42333 - 


■2 


1.40167 


1.41011 + 1 


1.87917 


1.67192 + 1 




Ae 


2.48012 - 


■2 


1.46078 


1.42938 + 1 


1.56431 


2.25171 + 1 




At 


1.88720 - 


■2 


1.50218 


1.44237 + 1 


1.35090 


2.92006 + 1 




As 


1.48800 - 


-2 


1.53215 


1.45153 + 1 


1.20040 


3.67684 + 1 




Ag 


1.20536 - 


-2 


1.55447 


1.45822 + 1 


1.09066 


4.52194 + 1 




Aio 


9.97399 - 


■3 


1.57151 


1.46325 + 1 


1.00834 


5.45528 + 1 


n = 2 


Ai 


7.45600 - 


-2 


1.73936 


2.70417 + 1 


1.54539 + 1 


7.51221 




A2 


4.92536 - 


-2 


1.96217 


2.77419 + 1 


1.48564 + 1 


1.06080 + 1 




A3 


3.28764 -^ 


-2 


2.12083 


2.82502 + 1 


1.43747 + 1 


1.49908 + 1 




A4 


2.39851 - 


-2 


2.22450 


2.85612 + 1 


1.40721 + 1 


1.98578 + 1 




As 


1.83252 - 


■2 


2.29718 


2.87742 + 1 


1.38614 + 1 


2.53951 + 1 




Ae 


1.44990 - 


-2 


2.35003 


2.89260 + 1 


1.37100 + 1 


3.15779 + 1 




At 


1.17785- 


-2 


2.38962 


2.90381 + 1 


1.35976 + 1 


3.84111 + 1 




As 


9.76979 - 


-3 


2.42001 


2.91234 + 1 


1.35120 + 1 


4.58946 + 1 




A9 


8.24157 - 


-3 


2.44384 


2.91897 + 1 


1.34452 + 1 


5.40287 + 1 




Aio 


7.05027 - 


-3 


2.46286 


2.92423 + 1 


1.33922 + 1 


6.28133 + 1 


n = 2.5 


Ai 


2.00471 - 


■2 


4.15142 


7.49546 + 1 


5.86265 + 1 


1.79915 + 1 




A2 


1.65192 - 


■2 


4.23203 


7.52209 + 1 


5.85001 + 1 


2.14747 + 1 




A3 


1.32031 - 


-2 


4.31214 


7.54815 + 1 


5.83730 + 1 


2.64277 + 1 




A4 


1.09699 - 


-2 


4.36892 


7.56633 + 1 


5.82826 + 1 


3.14401 + 1 




As 


9.23785 - 


■3 


4.41467 


7.58082 + 1 


5.82096 + 1 


3.69897 + 1 




Ae 


7.89994 - 


■3 


4.45110 


7.59226 + 1 


5.81513 + 1 


4.29372 + 1 




At 


6.83610 - 


-3 


4.48080 


7.60151 + 1 


5.81038 + 1 


4.93243 + 1 




As 


5.97686 - 


-3 


4.50527 


7.60909 + 1 


5.80646 + 1 


5.61402 + 1 




A9 


5.27200 - 


-3 


4.52569 


7.61538 + 1 


5.80320 + 1 


6.33882 + 1 




Aio 


4.68630 - 


-3 


4.54289 


7.62066 + 1 


5.80044 + 1 


7.10677 + 1 



Table 4. A table of integrals for the polytropic indices n - 
of the magnetic field. 



1, 1.5, 2 and 2.5. A^, which is given in Table 1, characterizes the configuration 
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